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FOREWORD 


This  report  describes  the  results  of  an  initial  investigation  of  new  methods  for  predicting  the 
future  trajectory  of  maneuvering  targets.  With  the  advent  of  highly  maneuverable  threats,  target 
trajectory  prediction  is  critical  for  the  success  of  ship  self-defense  systems.  This  report  discusses 
many  of  the  concepts  associated  with  trajectory  prediction  with  respect  to  the  dynamic 
constraints  of  the  target  and  its  self-imposed  limitations.  The  prediction  methods  have 
applications  in  command  and  decision  systems  and  weapons  control  systems,  such  as  the 
PHALANX  Close-In  Weapons  System,  NATO  Sea  Sparrow,  etc. 

Appreciation  is  given  to  Ernest  Ohlmeyer  of  the  Aeromechanics  Branch,  Naval  Surface 
Warfare  Center  Dahlgren  Division  (NSWCDD),  for  his  helpful  suggestion  to  use  proportional 
navigation  techniques  for  goal-based  trajection  prediction  for  maneuvering  targets. 

The  work  described  in  this  report  was  supported  through  the  Office  of  Naval  Research  and, 
more  specifically,  the  Surface  Weapons  Systems  Technology  Block  Program,  managed  at 
NSWCDD  by  Robin  Staton. 


Approved  by: 

D.  B.  COLBY,  Head 

Systems  Research  and  Technology  Department 
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ABSTRACT 


The  classical  problem  of  weapons  control  involves  predicting  the  future  position  of  a 
maneuvering  target.  Methods  of  target  trajectory  prediction  for  ship  self-defense  against  a 
maneuvering  antiship  cruise  missile  (ASCI^  are  studied.  Several  prediction  schemes  are 
described  and  studied  in  the  context  of  a  typical  single  ASCM  engagement.  Methods  of  utilizing 
the  estimated  threat  acceleration  and/or  a  "goal”  hypothesis  are  considered.  "Terminal 
Parameters"  (TP)  for  effectively  characterizing  the  threat  of  a  target  to  the  ship  are  introduced, 
alo^g  with  some  Measures  of  Effectiveness  (MOE)  for  comparing  the  different  prediction 
methods.  These  TPs  and  MOEs  are  illustrated  through  a  simulated  engagement  involving  a 
hypothesized  threat  track.  An  enhancement  of  defensive  capability  achieved  by  utilizing  a  goal 
hypothesis  is  demonstrated. 
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I 

I 

1 

CHAPTER  1 

'  INTRODUCTION 


One  of  the  most  important  and  difficult  problems  facing  the  surface  Navy  is  self-defense 
against  fast,  maneuvering  antiship  cruise  missiles  (ASCMs),  which  may  skim  the  sea  surface  aiwi 
have  a  small  radar  cross  section.  Among  the  many  tasks  that  need  to  be  completed  to 
successfully  address  this  problem  is  that  of  improving  techniques  for  threat-trajectory  prediction. 
The  ability  of  future  ASCMs  to  carry  out  sharp  and  deceptive  maneuvers  inteijects  much 
uncertainty  into  the  trajectory-prediction  process. 

This  report  describes  some  simple  algorithms  for  predicting  the  future  trajectory  of  an 
incoming  ASCM.  The  algorithms  are  separated  into  two  general  classes.  One  uses  only  the 
threat  state  estimates  (position,  velocity,  and  perhaps  acceleration)  at  the  instant  of  firing  and 
is  called  "State-Based. "  The  other  class  assumes  additional  knowledge  of  the  threat’s  destination 
or  goal  and  is  called  "Goal-Based. "  The  latter  method  is  more  pertinent  to  self-defense  than  to 
situations  involving  a  task  force  where  alternative  threat  destinations  are  feasible. 

The  threat  speed  for  each  of  the  new  predictors  is  assumed  constant  throughout  the  duration 
of  each  predicted  trajectory.  Several  effects  are  thereby  neglected:  energy  conversion  between 
potential  and  kinetic  caused  by  rising  or  descending,  the  usual  frictional  loss  of  energy,  and  the 
possibility  of  the  threat  varying  its  thrust  or  drag  coefficient  to  confuse  the  defenders.  While 
these  effects  can  be  important  in  some  circumstances,  it  is  worthwhile  to  gain  understanding  of 
the  simpler  constant-speed  case  before  introducing  complications.  For  a  threat  with  limited 
capability  of  controlling  its  longitudinal  acceleration,  and  whose  trajectory  is  not  very  steep, 
these  simplifications  are  appropriate.  In  the  present  state  of  ASCM  technology,  rapid  speed 
variation  is  not  feasible,  especially  in  the  endgame  portion  of  the  flight. 

State-Based  predictors  include  the  well-known  constant-velocity  (CV)  and  constant- 
acceleration  (CA)  predictors.  The  new  state-based  predictors  include  Constant  Turning  Rate 
(CTR),  Exponentially  Decreasing  Turning  Rate  (EDTR),  and  Helical.  The  CTR  predictor  is 
intended  as  a  replacement  of  CA  predictor,  which  does  not  maintain  constant  threat  speed.  The 
EDTR  predictor  transitions  from  CTR  to  CV,  and  thus  avoids  prolonged  continuation  of  a 
turning  maneuver  while  maintaining  constant  speed.  The  Helical  predictor,  which  produces  a 
constant-curvature  and  constant-torsion  trajectory,  is  proposed  for  cases  in  which  there  is  reliable 
threat  jerk  information.  The  trajectory  resulting  from  constant  lateral  and  longitudinal 
acceleration  (CLLA)  is  discussed.  In  this  trajectory  the  acceleration  vector,  loosely  speaking. 
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is  constant  when  referred  to  the  velocity  vector,  which  is  quite  different  from  the  CA  predictor, 
in  which  the  acceleration  vector  is  constant  in  an  inertial  frame. 

The  Goal-Based  predictors  include  the  Minimax  (MM),  the  Earliest- Arrival  (EA),  and 
several  Proportional-Navigation  (PN)  predictors.  The  MM  predictor  minimizes  the  maximum 
lateral-acceleration  magnitude  required  to  reach  the  ship.  TTie  EA  predictor  assumes  that  the 
threat  will  get  to  the  ship  as  soon  as  possible  while  maintaining  constant  speed.  The  PN 
predictors  steer  the  threat  toward  the  ship  through  PN  formulas  that  specify  an  appropriate 
acceleration  at  any  instant  (i.e.,  ignoring  the  initial-acceleration  information).  The  Twisted  PN 
(TPN)  predictor  is  proposed  for  cases  in  which  it  is  desirable  to  make  use  of  the  initial  estimated 
threat  acceleration.  The  TPN  predictor  produces  a  trajectory  with  nonvanishing  torsion. 

This  report  is  organized  as  follows.  The  analysis  uses  a  simple  vector  formalism  as 
introduced  in  Chapter  2.  Chapter  3  introduces  the  traditional  and  new  state-based  prediction 
methods.  The  goal-based  prediction  methods  are  described  in  Chapter  4.  Chapter  5  introduces 
the  concept  of  Terminal  Parameters  (TPs),  which  translate  the  initial  threat  state  estimates  into 
values  that  facilitate  the  command  and  control  process  (e.g.,  whether  to  fire,  when  to  fire,  how 
to  predict  the  threat  trajectory).  Measures  of  Effectiveness  (MOEs)  are  introduced  in  Chapter 
6,  along  with  a  comparison  of  the  methods  through  a  test  scenario.  Chapter  7  considers 
modifications  needed  when  the  predicted  trajectory  would  otherwise  take  the  threat  through  the 
sea  surface.  A  general  summary  of  the  work  and  concluding  remarks  are  given  in  Chapter  8. 
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CHAPTER! 

PROBLEM  FORMULATION 


PLANAR  TRAJECTORIES 

Let  the  threat  position,  velocity,  and  acceleration  at  time  t  =  0  be  denoted  by  , 

respectively.  Let  zit)  be  the  predicted  threat  position  based  on  these  values.  In  a  large 
number  of  cases,  it  is  desirable  to  assume  that  the  threat  will  continue  to  fly  in  the  plane 
osculating  to  the  true  trajectory  at  the  time  t  =  0 .  The  formula 

Z(t)  =  Xg  +  a  (OXjj  +  p  (2-1) 

ensures  that  the  threat  remain  on  this  plane.  In  this  case,  the  trajectory  is  completely  determined 
by  the  two  scalar  fimctions  a  ( t)  and  P  ( t)  .  Note  that  these  functions  and  their  derivatives 
are  chosen  to  have  the  mitial  values 

a(0)=0  p(0)=0 

d(0)=l  p(0)=0  (2-2) 

a(0)=0  P(0)=1 

For  constant-speed  motion,  the  constraint  relation 

=  \z,\^r  (2-3) 

is  required  by  (2-1). 


CONSTRAINT  ON  THE  THREAT  ACCELERATION 

For  the  prediction  schemes  discussed  in  this  report,  the  threat  speed  is  assumed  to  remain 
constant  throughout  the  prediction  time.  This  requires  that  the  acceleration  always  remain 
perpendicular  to  the  velocity.  If  the  initial  conditions  or  state  estimates  obey  this  constraint,  the 
prediction  formulas  satisfy  the  constraint  for  all  time.  When  the  initial  values  are  obtained  from 
sensor  information  by  means  of  some  track  filtering  scheme,  this  constraint  will  not  generally 
be  satisfied.  In  this  case,  one  of  the  two  vectors  can  be  altered.  Since  the  estimated  velocity 
value  is  more  reliable  than  the  acceleration  estimate,  the  acceleration  is  altered  when  necessary; 
assuming  that  the  estimated  acceleration  is  of  some  value,  it  should  be  altered  as  little  as 
possible.  The  following  criteria  were  considered  to  determine  a  modified  acceleration,  ,  to 
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replace  the  initial  estimated  acceleration. 


The  modified  acceleration  should  be  orthogonal  to  the  estimated  velocity.  The  modified 
acceleration  should  lie  in  the  plane  of  the  velocity  and  the  estimated  acceleration.  The 
magnitude  of  the  modified  acceleration  will  be  set  equal  to  the  component  of  the  estimated 
acceleration  that  is  perpendicular  to  the  velocity  and  given  by 


(2-4) 


where  is  the  estimated  velocity  vector  and  is  the  estimated  acceleration  vector  as  provided 
by  a  track-filtering  algorithm  such  as  a  Kalman  filter. 
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CHAPTERS 

STATE-BASED  PREDICTORS 


The  trajectory  predictors  can  conveniently  be  divided  into  two  classes;  state-based  and  goal- 
based.  The  state-based  predictors  depend  only  on  the  estimated  threat  state  (position,  velocity, 
and  perhaps  acceleration),  and  are  discussed  in  this  chapter,  while  the  goal-based  predictors  are 
discussed  in  Chapter  4. 


CV  TRAJECTORY 


The  CV  trajectory  is  given  by  (2-1),  where 


a  ( t)  =  t 
P(t)  =  0 


(3-1) 


The  CV  threat  continues  to  fly  in  a  straight  line  at  constant  speed.  The  constancy  of  the  speed 
is  demonstrated  by  the  fact  that  the  expressions  (3-1)  satisfy  the  condition  (2-3).  The  estimated 
initial  threat  acceleration  is  not  used  in  the  CV  predictor. 


CA  TRAJECTORY 


The  CA  trajectory  is  given  by  (2-1),  where 


a  (c)  =  t 
P(fc)  =  tV2 


(3-2) 


Here,  the  acceleration  vector  is  constant  in  magnitude  and  direction  if  referred  to  an  inertial 
frame.  Since  the  expressions  (3-2)  do  not  satisfy  the  condition  (2-3)  for  constant  speed,  the 
speed  generally  varies. 


CTR  TRAJECTORY 

In  order  to  achieve  a  trajectory  that  is  qualitatively  similar  to  CA,  but  which  obeys  the 
constant-speed  constraint,  the  CTR  predictor  is  used.  The  trajectory  is  given  by  (2-1)  with 

a(t)  =  a>o^sin(&)ot)  (3-3) 

P  ( t)  =  (Oq^  [1-cos  (Wot)  ] 
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where  Uq  is  the  constant  turning  rate  given  by 


(3-4) 


The  trajectory  is  of  course  a  circular  arc.  Since  the  expressions  (3-3)  satisfy  the  condition  (2-3), 
the  speed  remains  constant  over  the  time  of  prediction. 


EDTR  TRAJECTORY 

It  is  unreasonable  to  expect  that  a  nonzero  turning  rate  will  persist  indefinitely,  because  the 
threat  would  be  flying  in  circles.  The  EDTR  trajectory  allows  for  the  turning  rate  to  decay  to 
zero.  The  resulting  trajectory  begins  similar  to  a  circular  arc  and  ends  up  as  a  CV  trajectory. 
This  EDTR  trajectory  should  not  be  confused  with  the  linear  exponentially  decaying  acceleration 
predictor,  which  results  in  a  variable-speed  trajectory. 

The  constraint  (2-3)  is  satisfied  by  setting 

d  =  cos[p(t)],  0  =  Wo^sin  [p  ( t)  ]  (3-5) 

where  p  ( t)  is  an  undetermined  function.  No  essential  additional  assumption  has  been  made 
by  this  substitution.  Satisfying  the  initial  conditions  on  a  and  0  of  (2-2)  requires  that 

p(0)=0,  p(0)=a)o  (3-6) 

The  first  of  these  relations  is  obtained  by  setting  t  equal  to  zero  in  (3-5)  and  applying  the  initial 
conditions  (2-2).  The  second  is  obtained  by  setting  t  equal  to  zero  in  the  derivative  of  (3-5) 
and  applying  (2-2).  The  form  of  the  lateral  acceleration  needs  to  be  specified.  Differentiating 
the  position  vector  (2-1)  twice  and  taking  its  magnitude  gives 

!£|  =  |i|rol 

where  &  ( t)  and  0  ( t)  are  eliminated  by  using  (3-5).  For  the  exponentially  decreasing  turning 
rate,  we  take  an  exponentially-decreasing  lateral  acceleration  with  time  constant  t  and  set  the 
above  expression  equal  to 

\l\  =  (3-8) 

This  requires  that 

p  =  (dge'‘^/^  (3-9) 

This  differential  equation,  along  with  the  initial  values,  determines  p  to  be 

p  ( t)  =<i>oT  (3-10) 

The  fiiiKtions  a  and  0  can  now  be  obtained  by  integrating  their  previously  determined 
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derivatives  given  by 


a  ( t)  =  Jcos  [WqT  ]  dt 


(3-11) 


P(t)  =a>o^j*sin[<i)oT  ]  dt 


Making  the  substitution 


in  (3-11)  and  (3-12)  gives 


A,  =  6>oTe 


(3-12) 


(3-13) 


—  11  — 

a(t)=-T  J  cos  (o>oT-l) 


(3-14) 


—  u  - 

P(t)=-—  f  sinCWpT-A) 

(On  J 


(3-15) 


The  sine  and  cosine  integral  functions  can  be  used  to  express  the  solution  of  these  integrals.  In 
the  notation  of  Abramowitz  and  Stegun*  they  are  expressed  as 


S2(z)  =  f-^^dt=T 

J  t  _ n 


(_l)n^2n*l 


^0  (273+1)  (273+1)  ! 


(3-16) 


Ci  (z)  =  Y+lnz  +  J 


cost-1 


=  Y+ln^  *T 

^  273  (273)  ! 


(3-17) 


where  y  =  0 . 577 22  is  Euler’s  constant.  The  sums  in  the  expressions  for  Si  ( z)  and  Ci  ( z) 
converge  slowly  for  large  argument,  for  which  Abramowitz  and  Stegun  provide  good 
approximations  in  the  form  of  rational  algebraic  expressions.  The  result  can  be  expressed  in 
the  form 


a  ( t)  =  -tIcos  {(OqT)  [Ci  (tOgTe  -  Ci  ((OqT)  ] 

+  sin  (UqT)  [Si  i(0gze~'^^^)  -Si  (WqT)  ] } 


(3-18) 
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P{t)  =  TO>o^{ -sin  (Ogt)  [Ci  (o>oTe'^'^^)  -  Ci  (oJoT)  ] 

+  cos  ((OqT)  [Si  (OqT  -  Si  (Wq't)  ] } 


(3-19) 


Approximations  of  (3-18)  and  (3-19)  for  small  and  for  large  values  of  t  are  readily 
obtained.  For  t>T  (3-18)  and  (3-19)  approach  the  asymptotic  values 


and 


a  ( t)  -  -  xfcos  (WqT)  [y  +  In  (oJqT)  -  -  Ci  (tOoX)  ] 

+  sin  ((OqT)  [o)qT -  Si  (OqT)  ])  + 


(3-20) 


P  (t)  -  -TQoMsin  ((dot)  [y  +  ln  (WqT)  -  tf^-Ci  ((OqT)  ]  (3-21) 

-  cos  (OqT)  [<doTe''^/^-Si  (WoT)  ] }  + 

where  means  "order  of"  and  indicates  a  quantity  whose  ratio  to  remains 

bounded  as  &-<».  For  t<T  the  same  expressions  approach  the  asymptotic  values 


a  it)  =  t-  +  0(  t^)  , 

6 


P{t) -4---|^+0(£;^)  (3-22) 

2  6  T 


where  here  Oi  t^)  indicates  a  quantity  whose  ratio  to  remains  bounded  as  t-0 .  The  time 
derivatives  of  the  asymptotic  expressions  (3-20)  and  (3-21)  for  large  t  indicate  that  the  velocity 
vector  fmally  approaches 


=±oCos  (WqT)  +(i)o 

The  angle  between  the  initial  and  final  velocities  is  determined  by  the  relation 

t'tf 

-~.z  =  cos  (WoT) 


(3-23) 


(3-24) 


which  shows  that  the  total  change  in  heading  (amount  of  turning)  is  equal  to  UqT 


It  is  seen  that  for  positive  x  (decreasing  turning  rate),  the  velocity  vector  eventually 
approaches  the  final  constant  value  given  by  (3-23).  Thereafter,  the  trajectory  is  a  straight  line, 
which  can  be  represented  by 

x(  t)  =ro+2t+i:^t  (3-25) 

which  defines  the  constant  displacement  vector  h.  This  vector  is  evaluated  by  discarding  the 
exponentially  decreasing  terms  in  (3-20)  and  (3-21)  and  comparing  (2-1)  with  (3-25).  The 
result  is 
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h  =  tI-cos  ((OqT)  [y  +  ln  (Wo^)  -Ci  (cOqT)  ] 

+  sin  (OqT)  Si  (QoT)}i^ 

+  — {-sin  ((OqT)  [y  +  In  (OqT)  -Ci((i)oT)]  (3-26) 

-  cos  (  (OqT  )  Si  (  (OqT  )  IjEq 

The  length  of  this  vector  is 

\h\  =  Tlr^KCy  +  ln  (<i)(jT) -ci  (OqT)  ]  2  +  5i2  (0)qT)}^/2  (3-27) 

Figure  3-1  shows  both  the  actual  trajectory  (solid  curve)  given  by  (2-1)  and  a  hypothetical 
trajectory  (dashed  line)  given  by  (3-25)  representing  the  motion  of  a  hypothetical  missile  that 
travels  during  its  entire  flight  with  the  constant  velocity  equal  to  the  real  missile’s  final  velocity. 
This  hypothetical  missile  began  its  flight  from  the  initial  position  at  the  beginning  of  the  ray 
shown  on  Figure  3-1  at  the  same  time  the  real  missile  was  at  its  initial  position. 


X  (orbitrory  units) 


nCURE  3-1 .  TYPICAL  EDTR  TRAJECTORY  (SOUD)  AND  ASYMPTOTE  (DASHED)  AS 
RAY  DEPARTING  FROM  A  POINT  THAT  IS  SYNCHRONOUS  WITH 
INITIAL  THREAT  POSITION 
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Projected  along  its  linear  trajectory,  the  hypothetical  missile  begins  its  flight  a  distance  of 

b‘£ 

-r-r-j  =  -T|iJ  [y  +  ln  (OqT) -Ci  (WoX)  ]  (3-28) 

1^0 . 

behind  the  initial  position  of  the  real  missile  (but  eventually  catches  up  and  asymptotically 
coincides  in  position  with  the  real  missile).  The  "offset"  of  the  hypothetical  trajectory;  i.e.,  the 
minimum  distance  of  the  hypothetical  missile  from  the  initial  position  of  the  real  missile,  is 
given  by 

6  =  {|2i|2-  =  T|rJSi(o)oT)  (3-29) 


An  appropriate  value  of  the  die-off  time  constant,  t  ,  must  be  selected.  If  this  value  greatly 
exceeds  die  prediction  time,  the  EDTR  will  not  perceptibly  differ  from  CTR.  It  may  be  easiest 
to  base  the  value  of  x  on  the  expected  total  change  of  heading,  (OqT  . 


HEUCAL  TRAJECTORY 

The  helix,  being  a  curve  having  constant  curvature  and  constant  torsion,  is  the  simplest 
reasonable  nonplanar  trajectory.  It  requires  loiowledge  of  the  threat  jerk,  in  addition  to  the 
position,  velocity  and  acceleration,  from  the  estimated  state  vector.  Although  the  jerk  estimate 
may  be  rather  poor,  it  is  constrained  to  improve  prediction  in  a  manner  to  be  described  below. 
It  should  be  noted  that  prediction  models  based  on  planar  motion  tacitly  assume  that  the  jerk  lies 
in  the  velocity-acceleration  plane.  This  assumption  may  be  worse  in  some  instances  than  using 
jerk  information  from  the  estimated  state  vector. 


The  helix  is  described  by  a  simple,  well-known  formula  that  expresses  the  position  vector 
as  a  function  of  time.  The  constants  of  this  representation  are  expressed  as  the  initial  position, 
Xg,  velocity,  x^,  acceleration,  and  jerk,  Xg.  The  helical  trajectory  is  represented  by 


where 


X(t)  =  Xg+  {Xg  +  «0^Xg)  t 

+  &)o^  [1-COS  (Wgt)  ]Xg-(*>o^XgSin  (Wgt) 


(3-30) 


Since  for  constant  speed  the  acceleration  is  always  perpendicular  to  the  velocity. 


(3-31) 
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if  =  0  (3-32) 

Taking  the  derivative  of  (3-32),  it  is  seen  that 

=  0  (3-33) 

which  indicates  that  the  angle  between  the  jerk  and  the  velocity  lies  between  71/2  and  371/2 . 
Taking  derivatives  of  (3-30)  for  direct  substitution  into  (3-32),  using  (3-31)  and  (3-33),  it  is  seen 
that 


i-f  =  [cos  (Wot) -cos  {2&>ot)  ]  «0  (3-34) 

This  shows  that  the  initial  jerk  must  be  perpeiKiicular  to  the  initial  acceleration  in  order  that  the 
expression  (3-34)  vanish  for  all  time.  In  a  similar  manner  it  is  found  that 

f-x  =  lo'-^oCOS  (2a>ot)  (3-35) 

This  shows  that  the  jerk  is  always  perpendicular  to  the  acceleration,  since  initially  this  is  so. 
Since  the  left  side  of  (3-35)  is  (halO  derivative  of  the  acceleration  dotted  onto  itself,  the 
acceleration  magnitude  remains  constant.  Also,  since  Wq  remains  constant,  (3-31)  requires  that 
the  jerk  magnitude  remain  constant. 

Note  that  the  helical  trajectory  reduces  to  the  familiar  constant  turning  rate  (in  a  circular 
trajectory)  if 


(3-36) 


This  suggests  that  it  may  be  convenient  to  represent  the  jerk  as  a  deviation  from  that  required 
to  sustain  circular  motion.  Accordingly,  the  vector  p  is  defmed  as  this  deviation  and  we  have 

p(t)  =  f  (t)  +<i>of  ( t)  (3-37) 

Thus,  (3-36)  and  (3-37)  show  that  p(t)  -0  in  the  special  case  where  the  helical  motion 
degenerates  to  pure  circular  motion.  The  right  side  of  (3-37)  can  be  evaluated  by  taking  time 
derivatives  of  (3-30)  to  show  that  the  vector,  p,  remains  constant.  It  is  also  seen  that  p  is 
parallel  to  the  nonfluctuating  part  of  £■  Thus,  p  is  parallel  to  the  axis  of  the  helix. 
Accordingly,  we  write 

p(t)=p„  (3-38) 

The  position  vector  can  then  be  written  as 

r(t)  =  Xo  +  «o^Pot+ («o^Xo-«o^Po)sin(&)ot)  (3-39) 

+  0)o^Il-COS((i>ot)]f^ 


Also, 
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i-Po  =  “o^IPqP  (3-40) 

follows  from  (3-37),  (3-38),  and  (3-39),  showing  that  the  angle  between  the  velocity  vector  and 

As  noted,  the  initial  estimated  values  of  the  threat’s  state  vector  must  satisfy  certain 
constraints.  As  in  the  case  of  the  threat  acceleration  vector,  the  jerk  must  also  be  made  to 
satisfy  some  constraints.  If  no  meaningful  estimate  of  the  jerk  is  available  then  there  is  no 
reason  to  use  the  helical  prediction.  Accordingly,  it  is  asstimed  here  that  reasonable  values  of 
the  threat  position,  velocity,  acceleration,  and  jerk  are  available.  As  the  estimated  acceleration 
value  is  more  reliable  than  the  jerk,  the  acceleration  estimate  will  be  adjusted  first  by  replacing 
its  estimated  value  by  defmed  in  (2-4).  To  find  a  suitable  replacement,  for  the 
estimated  jerk  vector,  the  following  criteria  are  considered. 

(a)  The  resulting  jerk  vector  should  be  perpeiKiicular  to  the  acceleration;  i.e. ,  =  0 . 

(b)  The  relation  =  0,  as  required  by  (3-33). 

(c)  The  should  be  equal  to  the  component  of  which  is  perpendicular  to  . 

The  first  two  criteria  are  essential.  It  turns  out  that  in  general  the  third  cannot  be  satisfied 
exactly.  It  may  be  desirable  to  satisfy  criterion  (c)  as  well  as  possible,  in  a  least-square-error 
sense. 


Let  us  represent  the  vector,  which  is  to  replace  the  estimated  jerk,  by 

which  satisfies  criterion  (a)  for  all  values  of  p  and  v ,  which  are  scalar  quantities  to  be 
determined.  Criterion  (b)  is  satisfied  by  requiring 


M  = 


(3-42) 


The  value  of  v  is  still  to  be  determined.  It  is  noted  that  so  far  no  use  at  all  has  been  made  of 
the  estimated  jerk  vector.  Criterion  (c)  is  relaxed  to  minimize  (with  respect  to  v)  the 
discrepancy  between  x^,^  and  the  component  J?  of  x^  which  is  perpendicular  to  This 
component  is  equal  to 


Y  = 


X  ’£ 

s  _  0  OiH  a 

*0  la  12 


(3-43) 
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Thus,  "equal"  is  replaced  by  "close  as  possible"  in  a  least-square-error  sense  whereby  the 
quantity  Yj  is  minimized  with  respect  to  v  .  This  determines  v  to  have  the  value 


V 


and  the  modified  jerk  is  now  determined  by  (3-41). 


(3-44) 


CONSTANT  LATERAL  AND  LONGITUDINAL  ACCELERATIONS 


Let  aj  be  the  magnitude  of  the  constant  longitudinal  or  thrusting  acceleration,  and  let 
be  the  magnitude  of  the  constant  lateral  acceleration.  Only  the  magnitudes  are  constant.  The 
vector  accelerations  will  be  denoted  by  a ^  and  respectively,  and  are  not  constant.  Gravity 
is  neglected,  and  the  trajectory  is  assumed  to  lie  in  a  fixed  plane,  the  "flight  plane" .  Note  that 
if  the  flight  plane  is  horizontal,  this  assumption  is  not  needed.  Let  this  plane  be  determined  by 
a  normal  unit  vector,  n.  Let  z(t)  be  the  position  at  time  t.  According  to  these  definitions. 


^ 


(3-45) 


and 


where 


A  ’A 


=  0 


The  longitudinal  acceleration  vector  is 


w 


Let  the  lateral  acceleration  vector  be  represented  by 


(3-46) 


(3-47) 


A^  =  \iaxi  (3-48) 

where  p  is  a  scalar  to  be  determined.  This  form  ensures  that  the  lateral  acceleration  will  lie  in 
the  flight  plane  and  be  perpendicular  to  the  velocity  vector.  Dotting  this  expression  onto  itself 
gives 


=  af  =  p2[j^|2_  {0-^)2]  (3-49) 

which  allows  p  to  be  expressed  in  terms  of  other  quantities.  The  last  term  within  the  brackets 
on  the  right  is  zero,  because  the  velocity  vector  lies  in  the  flight  plane.  Substituting  this  value 
of  p  back  into  the  expression  for  the  lateral  acceleration  (3-48)  gives 


A 


L 


inxi) 

i 


(3-50) 
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The  total  acceleration  is  the  sum  of  the  longitudinal  and  lateral  accelerations,  in  accordance 
with  (3-46),  and  we  have 


£  = 


a^(nx£)  +aj£ 


(3-51) 


Let  a  cartesian  coordinate  system  on  the  flight  plane  be  established  with  unit  vectors  £  and  £ 
mutually  perpendicular  and  perpendicular  to  a.  Let  u  and  v  be  the  velocity  components  in  this 
system,  according  to 

£=ul*vi  (3-52) 


In  accordance  with  (3-51),  the  acceleration  components  are  given  by 


and 


-a.v+a-u 
u  =  — ^ 


(3-53) 


Equations  (3-53)  and  (3-54)  comprise  a  system  of  first  order  differential  equations  in  the  velocity 
components.  To  find  a  solution,  it  is  convenient  to  rewrite  them  in  terms  of  the  speed,  w,  and 
direction  0.  Accordingly, 


U  =  WCOS0 
V  =  wsin0 

Substituting  these  expressions  into  (3-53)  and  (3-54),  we  obtain 

IVCOS0- w0sin0  =  -a^sin0  +  aj,cos0 
wsin0  +  w0cos0  =  a^cos0  +  aj.sin0 

Multiplying  these  equations  by  cos0  and  sin0,  respectively,  and  adding,  we  obtain 

w  =  aj.  (3-57) 

Similarly,  multiplying  the  two  equations  of  (3-56)  successively  by  -sin0  and  cos0  and 
adding,  we  obtain 

w0  =  a^  (3-58) 

Equations  (3-57)  and  (3-58)  show  that  the  simple  effects  of  the  two  acceleration  components  are 
more  easily  seen  in  terms  of  the  speed  and  direction  than  on  the  rectilinear  velocity  components. 
The  integrals  of  (3-57)  and  (3-58)  are 

w  =  Wg  +  aj-t  (3-59) 


(3-55) 

(3-56) 


and 
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J=  Jd+AtO  sin[5ln(l+AtO  ]  dt^ 


(3-68) 


where 


A  = 


)fuQ 


B  = 


2  2 


A  simple  substitution. 


A.=  — Ind+AtO 


allows  straightforward  evaluation  of  the  integrals.  The  results  are 

I  j  I  -Sr  jj 


^  ^  ^Uq  -t- Vq  I 


^  2  2 
4ar+a£ 


2ajjcos 


and 


The  coordinates  can  then  be  expressed  as 


and 


(t)  =  judt  =  u^I-v^J+Xo 


y(t)  =  fvdt  =  v^I+u^J+yo 


(3-69) 


(3-70) 


(3-71) 


(3-72) 


(3-73) 


(3-74) 


It  is  convenient  to  refer  the  cartesian  components  to  a  system  that  is  determined  by  the  initial 


and  .  Let  the  unit  vectors  1  and  i  be  defined  by 

^0  0 


1- 

W' 


^  ■  liJ 


(3-75) 


where  a  is  the  unit  vector  perpendicular  to  both  and  ,  given  by 


i  2  f  P  -  (i  '1  ) 


(3-76) 


The  acceleration  components  a^.  and  a^  are  also  obtained  from  the  initial  velocity  and 
acceleration  by  taking 
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and 


aj,  = 


I  t 

0  *-0 


(3-77) 


a 


L 


(3-78) 


The  above  formulas  for  a^  and  a^  simply  perpetuate  the  values  they  apparently  had  at  time  . 


The  prediction  algorithm  can  be  summarized  as  follows.  For  a  given  initial  state,  compute 
a^  aiKi  a^  from  (3-77)  and  (3-78),  the  unit  vector  a  from  (3-76),  and  the  unit  cartesian  vectors 
^  and  i  from  (3-75).  Then,  for  any  desired  time  t  the  position  vector  is  given  by 
x(  t)  =x{t)  ±  +  y{t) i,  where  x{  t)  and  y(  t)  are  given  by  (3-73)  and  (3-74)  with  I  aiKi 
J  given  by  (3-71)  and  (3-72), 
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CHAPTER  4 

GO  '4L-ORIENTED  APPROACH  TO  PREDICTION  FOR  SELF-DEFENSE 


The  goal-based  approach  predicts  future  threat  behavior  on  the  assumption  that  the  threat  is 
targeted  on  the  ship.  If  this  assumption  is  not  valid,  the  ship  is  safe  although  some  defensive 
rounds  may  be  wasted.  Clearly,  for  self-defense,  failure  to  make  this  assumption  could  be  risky. 

Since  not  all  threat  states  represent  a  danger  to  the  ship,  the  fust  section  of  this  chapter, 
"Domain  of  Threat  Capability,"  attempts  to  exclude  goal-based  prediction  from  those  targets 
whose  state  vector  indicates  that  the  ship  is  an  unlikely  target.  The  next  section,  "Constraint  on 
Initial  Acceleration, "  considers  the  issue  of  constraining  the  acceleration  estimate  so  that  the 
target  with  a  maximum  acceleration  limit  will  intercept  the  ship.  Because  the  threat  acceleration 
can  exhibit  quick  changes,  it  is  safer  to  assume  a  value  that  will  drive  the  threat  toward  the  ship. 
The  remainder  of  the  chapter  discusses  several  particular  models  for  prediction,  each  of  which 
embodies  the  goal-oriented  feature  through  different  criteria. 


DOMAIN  OF  THREAT  CAPABIUTY 

It  is  reasonable  to  assume  that  a  threat  in  a  given  initial  state  would  be  capable  of  reaching 
any  ship  within  a  certain  region  and  unable  to  reach  a  ship  outside  this  region.  Inverting  this 
concept,  for  a  given  ship  position,  there  are  some  target  states  that  would  be  capable  of  hitting 
the  ship,  while  all  other  states  would  not.  For  a  given  initial  target  state  and  ship  position,  it 
is  iKcessary  to  determine  fust  if  the  threat  poses  any  danger  to  the  ship.  Excluding  horn 
consideration  those  states  considered  not  dangerous  implies  some  limitations  on  threat 
maneuverability.  One  such  limitation  might  be  to  assume  that  it  cannot  overfly  the  ship  and 
return  to  hit  it.  The  following  restrictions  attempt  to  formalize  such  limitations  in  a  simple  way. 
The- threat  is  assumed  unable  to  reach  any  point  in  the  halfspace  behind  itself.  It  is  also  assum^ 
that  the  theat  is  unable  to  reach  any  point  within  the  "exclusion"  torus,  which  is  centered  on  the 
threat’s  initial  position  and  has  its  axis  of  symmetry  coincident  with  the  velocity  vector.  The 
radius  of  die  torus  generating  circle  is  equal  to  the  radius  of  curvature  defined  by  the 
hypothesized  maximum  lateral  acceleration,  .  That  is. 


Furthermore,  the  distance  from  the  axis  of  symmetry  to  the  center  of  the  generating  circle  is  also 
jbj .  This  is  the  limiting  case  of  a  torus  on  the  verge  of  intersecting  itself.  All  of  the  space  not 
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in  the  threat’s  rear  halfspace  nor  in  its  exclusion  torus  is  considered  to  comprise  its  domain  of 
capability.  The  sign  of  determines  if  the  ship  is  in  the  excluded  halfplane.  To  determine 
whether  or  not  the  ship  lies  within  the  exclusion  torus,  the  minimax  acceleration  a„,  to  be 
described  below,  is  computed,  and  if  a^>  the  ship  is  within  this  torus. 


CONSTRAINT  ON  INITIAL  ACCELERATION 

For  the  goal-oriented  hypothesis,  the  initial  threat  acceleration  is  constrained  so  that  the 
trajectory  plane  includes  the  origin.  This  is  quite  different  from  the  acceleration  constraint 
discussed  in  Chapter  2.  Let  the  adjusted  acceleration  vector,  to  be  substituted  for  the  estimated 
state  acceleration,  be  denoted  by  This  new  acceleration  vector  satisfies  the  following 

criteria: 

1.  The  adjusted  acceleration  is  perpendicular  to  the  velocity.  If  represents  the  adjusted 
acceleration,  then 

=  0  <4-2) 

2.  The  adjusted  acceleration  lies  in  the  plane  of  the  threat  position  and  velocity.  This 
condition  is  applied  by  setting 

<4-3) 

where  p  and  v  are  scalar  constants  to  be  determined. 

3.  The  magnitude  of  the  adjusted  acceleration  is  equal  to  the  magnitude  of  the  estimated 
acceleration  component  that  is  perpendicular  to  the  velocity.  That  is, 

=  1^0- 1-^0  <4-4) 

The  vector  within  the  outer  bars  on  the  right  side  is  the  r  3mponent  of  the  estimated  acceleration 
perpendicular  to  the  velocity. 

These  three  conditions  are  sufficient  to  determine  the  values  of  p  and  v  except  for  sign 
ambiguity  which  is  resolved  by  requiring  that  the  adjusted  acceleration  be  directed  toward  the 
origin,  or  ^  0 .  Dotting  to  both  sides  of  (4-3)  and  applying  (4-2)  allows  one  of  the 

two  constants  to  be  expressed  in  terms  of  the  other.  For  example,  we  can  set 
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and  express  by 


*.s  “  t* 


I*  12  '*■0 


Squaring  (4-4)  gives 
Equating  this  with  the  square  of  (4-6)  gives 


\l  P  = 


1  P  -  i:  -^(±  ■£  )^ 
^0 '  ^0  '^^0  0 


\l  =  ± 


li*„P-  (r -xj 


1/2 


■■0  '  '“‘0  '  “0  ““O 

Substituting  this  value  into  (4-6),  with  the  appropriate  sign  of  p ,  gives 


1 

2 

-xl 

|X„P|X,P-(X„*X„)^ 

(4-6) 


(4-7) 


(4-8) 


(4-9) 


The  magnitude  of  this  modified  acceleration  is  consistent  with  (4.7)  and  never  exceeds  |Xg  | . 
The  sign  of  (4-9)  was  chosen  so  that  the  resulting  trajectory  will  be  concave  toward  the  origin. 


MM  PREDICTOR 


"Minimax"  indicates  the  predicted  threat  track  that  minimizes  the  maximum  lateral 
acceleration  needed  to  reach  the  origin  after  passing  through  a  terminal  CV  segment  of  specified 
length,  s.  It  is  not  surprising  that  minimax  requires  that  the  initial  trajectory  segment  be  a 
circular  arc.  The  entire  track  lies  in  the  plane  (i.e.,  the  plane  of  the  estimated  initial 

threat  position  and  velocity).  If  the  initial  distance  |Xq|  is  less  than  s,  then  in  general,  no 
minimax  track  is  possible. 


It  is  convenient  to  define  the  position  of  the  center  of  curvature  of  the  circular-arc  segment 
at  time  by  the  expression 


At  this  time  the  radius  vector  is 


PXg-  (p-1) 


(4-10) 
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X  -X 


=  (ji-1) 


I  A  12^0 


(4-11) 


which  lies  in  the  (x^ ,  )  plane  and  is  seen  to  be  perpendicular  to  x^ .  Here,  p  is  a  constant 

scalar  to  be  determined.  The  condition 


|x^|2  =  (4-12) 

is  a  statement  of  the  Pythagorean  theorem  in  the  right  triangle  defined  by  the  origin,  the  point 
where  the  circular  arc  transitions  to  the  CV'  segment,  and  the  position  of  x^.  Substituting  (4-10) 
into  (4-12)  gives  a  linear  equation  in  p  (because  the  quadratic  terms  cancel)  from  which  the 
value 


|xJ^(s"-»-|xJ^)  -2(Xp-X„)" 
2[|xJ"|ioP-{X,-Xo)^] 


(4-13) 


is  obtained.  To  determine  interception  parameters  (PIP,  time  of  flight,  etc.),  it  is  convenient 
to  frrst  calculate  6 ,  the  arc  length  traversed  during  the  circular  flight  segment.  If  a  is  the  angle 
at  Xp  subtended  by  the  CV  flight  segment,  then  a  +  0  is  the  angle  at  the  threat  between  its 
velocity  vector  and  the  origin.  The  CV  segment  and  x^  form  the  two  legs  of  a  right  triangle, 
from  which  it  is  seen  that 


There  is  also  the  relation 


cos  a 


cos  (a+0) 


(4-14) 


(4-15) 


Together  (4-14)  and  (4-15)  serve  to  determine  0.  The  turning  rate  within  the  circular  arc,  o) , 
is  given  by  (3-4).  The  time  spent  in  the  circular-arc  segment  is  0/o) ,  which  thus  occurs  during 
the  interval 


to  t  ^  to  + 


b> 


(4-16) 


EA  PREDICTOR 

It  is  assumed  that  the  threat  will  immediately  achieve  a  lateral  acceleration  in  the  (x^  /  x^ ) 
plane,  with  maximum  allowable  magnitude  being  the  "limiting  acceleration"  ,  which  must  be 
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Specified,  and  directed  to  approach  the  ship.  The  predicted  track  will  follow  the  resulting 
circular  arc  until  the  velocity  vector  points  toward  the  origin,  whereupon  a  terminal  CV  track 
is  taken.  The  formalism  of  the  MM  predictor,  discussed  above,  is  applicable  except  that  a 
different  p  is  needed.  Equating  the  two  alternative  expressions  for  the  radius  of  curvature  gives 


Then  substituting  from  (4-10)  gives 


(4-17) 


=  1- 


(4-«> 

instead  of  (4-13).  The  simple  track  just  described  gives  the  EA  predictor,  provided 


which  implies  that 


^  \r. 


(4-19) 


sin^e  ^  (4-20) 

2  (1-p) 

where  e  is  the  angle  between  and  This  relation  is  obtained  by  eliminating  in 
terms  of  e  in  (4-18),  substituting  (4-10)  and  (4-17)  into  the  inequality  (4-19),  and  eliminating 
IjCq  I  between  the  two  resulting  relations.  An  interception  is  still  possible  if  this  condition  is 
violated,  but  this  requires  a  more  complicated  track,  and  for  the  moment  will  not  be  considered 
further. 


PN  METHODS 

Proportional  Navigation  (PN)  is  often  used  to  guide  a  pursuing  object  (a  missile)  to  intercept 
a  pursued  object  (the  target).  The  required  missile  acceleration  is  expressed  as  a  function  of  the 
missile  and  target  positions  and  velocities.  Let  and  x^  denote  the  missile  and  target 
positions  in  the  usual  PN  formulation,  with  the  missile-to-target  relative  position  vector  being 
X=Xj.-x^.  Here  the  target  position  is  replaced  by  the  ship  position  at  the  origin  and  the  missile 
becomes  the  threat.  Thus,  x^  is  replaced  by  fl,  and  x^  is  replaced  by  -x. 

One  particular  PN  law,  designated  here  as  Standard  PN  (SPN),  has  been  the  subject  of  a 
number  of  studies.  It  considers  the  rotation  rate  of  the  Line  of  Sight  (LOS)  and  specifies  an 
acceleration  that  will  asymptotically  reduce  this  rotation  rate  to  zero.  A  parameter  n,  called  the 
PN  gain,  determines  the  speed  at  which  the  LOS  rotation  rate  is  nulled.  The  acceleration  is 
taken  perpendicular  to  the  LOS,  and  not  necessarily  perpendicular  to  the  velocity.  Accordingly, 
the  speed  in  general  is  not  constant.  This  is  referred  to  as  "True  PN, "  in  contrast  with  "Pure 
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PN,"  which  conserves  missile  speed. ^  The  formula  for  the  control  acceleration  for  pure  PN  is 
given  by 


(4-2i) 

A  Pure  PN  can  be  obtained  from  True  PN  by  taking  as  the  required  acceleration  just  the  c 
component  of  SPN  which  is  perpendicular  to  the  velocity.  Here,  MPN  represents  this 
modification  of  SPN.  The  formula  for  the  control  acceleration  for  MPN  is  given  by 

[Wr-  (r-r)x]  (4-22) 

Adler^  proposed  a  convenient  three-dimensional  (Pure)  PN  law,  which  here  is  designated  as 
APN  and  in  our  notation  is  written 

(r-r)x]  (4-23) 

where  n  is  again  a  numerical  gain,  which  must  be  specified.  Adler  showed  that  a  value  ofn 
greater  than  two  is  required  for  practical  intercept  trajectories.  The  value  n  =  4  gives 
reasonable  trajectories  for  both  MPN  and  APN.  As  both  PN  formulas  determine  a  required 
acceleration  in  the  plane  of  the  threat  velocity  and  position,  the  resulting  trajectory  lies  in  a  plane 
that  includes  the  origin. 

The  PN  lateral  acceleration  magnitude  monotonically  decreases  if  the  pursued  object  is 
stationary,  as  the  ship  in  the  situation  under  discussion.  Therefore,  if  the  threat  is  initially 
capable  of  achieving  the  required  PN  acceleration,  it  is  reasonable  to  assume  that  the  required 
acceleration  remains  achievable.  When  the  acceleration  is  initially  unachievable,  two  obvious 
strategies  suggest  themselves.  One  is  for  the  threat  to  fly  a  CTR  segment  at  maximum  allowed 
lateral  acceleration  until  the  PN  acceleration  becomes  achievable.  Another  is  to  adjust  n 
downward.  In  the  latter  case,  an  appropriate  value  of  n  is  limited  by 

n  <:  [|xP|x|^- (X'X)^] (4-24) 

for  Adler  PN. 

As  indicated  in  Table  4-1,  some  closed-form  solutions  to  PN  trajectories  and  range-time 
histories  have  been  obtained,  in  the  special  case  of  a  stationary  target. 


TWISTED  PN  (TPN)  PREDICTOR 

The  TPN  trajectory  is  not  necessarily  confmed  to  a  single  plane.  If,  at  the  moment  of 
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TABLE  4-1.  EXISTENCE  OF  PN  SOLUTION 


True  or 

Pure  PN 

Trajectory 

Solution 

Range  vs 

Time 

SPN 

True 

a 

a 

MPN 

Pure 

d 

b 

APN 

Pure 

d 

c 

TPN 

Pure 

b 

b 

Notes: 

(a)  Cochran  et  al^  give  solutions  for  n=3  and  n=4  for  general  CV  target  motion. 

(b)  No  known  solutions  are  available. 

(c)  Groves  and  Gray^  give  solutions  for  n=3  and  n=4 

(d)  Groves  and  Gray*  give  solutions  for  general  n. 


observation,  the  threat  is  moving  in  the  plane  of  its  estimated  velocity  and  acceleration,  it  cannot 
simultaneously  remain  on  this  plane  and  also  intercept  the  ship  if  the  ship  does  not  lie  on  the 
plane.  The  TPN  predictor  gradually  rotates  the  osculating  plane  of  threat  motion  until  it  contains 
the  target  ship. 

The  common  PN  predictors  determine  a  required  acceleration  f  based  on  current  values 
of  the  estimated  threat  position  and  velocity.  They  assume  that  the  acceleration  can  be  abruptly 
changed  at  will.  The  resulting  trajectory,  of  course,  is  torsionless  (planar).  Here,  a  slight 
generalization  of  the  usual  PN  procedure  is  used.  It  will  he  assumed  that  the  jerk,  but  not  the 
acceleration,  can  be  assigned  at  will. 

Any  suitable  PN  trajectory  diat  is  targeted  on  the  origin  can  be  used.  However,  instead  of 
using  immediately,  the  cur’-ent  acceleration  is  to  be  ‘moved’  gradually  to  the  desired  value 
One  of  the  simplest  ways  of  achieving  such  a  gradual  change  in  the  acceleration  is  to  apply 
a  jerk  that  is  proportional  to  the  difference  between  the  current  and  the  desired  acceleration. 
The  applied  jerk  must  satisfy  the  basic  constraint 

i:-x+|l|2  =  0  (4-25) 

(obtained  by  differentiating  lip  =  twice).  A  convenient  prediction  formula  is  obtained  by 
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taking 


X  =  T'Mf  -1)--^!^  (4-26) 

*•  |j^|2 

This  expression  satisfies  the  basic  constraint  (4-25)  and  tends  to  null  the  acceleration  discrepancy 
as  time  elapses.  The  parameter  t  is  a  time  constant,  which  can  be  assigned  somewhat 
arbitrarily,  but  controls  the  rate  at  which  the  acceleration  reaches  the  desired  level.  When  the 
acceleration  becomes  essentially  equal  to  that  required  by  PN  the  jerk  becomes 


X  =  - 


(4-27) 


and  is  directed  in  the  opposite  direction  of  the  velocity.  The  motion  will  henceforth  be  planar. 
Momentarily  the  threat  is  in  a  circular-arc  trajectory,  as  can  be  seen  by  the  following  argument. 
A  solution  of  (4-27)  for  which  both  |£|  and  |x|  are  constants  is 


I  = 


(X,-X) 


(4-28) 


where  x^  is  a  constant.  It  is  seen  that  (4-28)  represents  an  acceleration  directed  always  toward 
the  point  x^,  and  has  magnitude  consistent  with  circular  motion,  and  that  |1|  and  |x|  are 
constant  as  required.  Of  course,  this  motion  is  hypothetical  because  the  resulting  acceleration 
would  continu^ly  stray  from,  and  need  to  be  brought  back  toward  that  required  by  PN. 


The  jerk  formula  (4-26)  contains  two  parameters  {n  and  t  )  that  need  to  be  specified.  If 
the  time  constant  t  is  chosen  too  large,  the  resulting  predicted  trajectory  will  miss  the  origin. 
It  seems  reasonable  to  choose  t  to  be  some  fraction  of  the  remaining  time  to  intercept.  As  this 
time  continually  decreases,  it  is  worthwhile  to  use  a  time-varying  t .  One  convenient  possibility 
is  to  take 


T  (  t) 


kl 

m\t\ 


(4-29) 


where  m  is  a  selectable  nondimensional  parameter  that  determines  the  rate  at  which  the  required 
acceleration  approaches  the  value  given  by  the  PN  formula.  This  ensures  that  the  trajectory  pass 
through  the  origin.  The  TPN  formula  (4-26)  then  becomes 


^  mn\i\^ 

Am  i  .  A 


inn\t\ri  111 


ni|i:|  _ 
\£  --r^i 


(4-30) 


ixr  \r\ 

for  APN,  while  for  MPN,  simple  obvious  modifications  need  to  be  made.  However,  the  two 
parameters,  n  and  m,  must  be  specified. 
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CHAPTER  5 

TERMINAL  PARAMETERS 


Decisions  and  hypotheses  are  to  be  based  on  the  latest  available  values  of 
However,  it  is  convenient  to  translate  these  into  other  parameters  more  closely  related  to  the 
decision-making  process.  These  subsequent  parameters  are  used  to  facilitate  conceptualizing  the 
possible  future  threat  trajectories.  It  is  emphasized  that  these  Terminal  Parameters  (TPs)  depend 
only  on  the  threat  state. 

Some  of  the  TPs  depend  on  an  assumed  lateral  acceleration  magnitude  that  the  threat  is 
incapable  of  exceeding.  This  value  may  be  based  on  the  identity  of  the  threat,  or  may  be 
hypothesized.  The  larger  the  value  of  aj,  the  more  difficult  the  prediction  problem  becomes. 

The  TPs  are  based  on  the  supposition  that  the  initial  guess  of  the  threat  acceleration  is 
so  poorly  determined  that  the  hypothesis  1^=0.  is  as  good  as  any  other.  The  proposed  TPs 
defined  in  the  following  list  will  be  considered. 

-  The  MM  acceleration 

-  The  MM  arrival  time 

tgj  ~  The  earliest  arrival  time 

Yac,  "  Look  angle  at  acquisition 

—  The  latest  arrival  time 

The  MM  acceleration  is  that  required  to  achieve  a  MM  trajectory,  described  in  Chapter  4. 
Note  that  its  value  depends  on  the  parameter  s,  the  length  of  the  final  CV  flight  segment.  The 
t„-  is  the  total  flight  time  in  the  CTR  and  the  CV  segments.  The  is  the  total  flight  time 
assuming  that  the  threat  will  fly  the  EA  trajectory,  as  described  in  Chapter  4. 

The  definition  of  the  look  angle  at  acquisition,  Yacg>  requires  a  hypothetical  defensive 
missile  or  bullet  to  be  launched  against  the  tlu’eat.  "Acquisition"  then  means  "of  the  threat  by 
the  defensive  missile,"  and  the  "look  angle"  is  between  the  LOS  from  defensive  missile  to  threat 
and  the  velocity  of  the  defensive  missile.  A  small  look  angle  is  of  course  advantageous  to  the 
defensive  side.  In  the  examples  presented  here,  the  defensive  missile  is  a  bullet  assumed  to  fly 
a  CV  trajectory.  This  greatly  simplifies  the  computations.  However,  the  same  concepts  exist 
for  defensive  missiles  flying  other  types  of  trajectories.  While  the  criterion  for  acquisition  here 
is  artificially  simple,  it  could  have  been  made  to  depend  on  a  large  number  of  environmental 
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factors  and  technical  details  of  the  defensive  missile. 

The  latest  arrival  (LA)  time  is  useful  in  some  contexts,  but  its  definition  is  not  as 
circumscribed  as  the  EA  time.  For  this  reason,  tl^re  was  no  discussion  of  an  LA  trajectory  in 
Chapter  4.  The  LA  time  is  based  on  the  threat  continuing  a  reasonable  trajectory  that  delays  as 
long  as  reasonably  practical  before  turning  toward  the  ship.  This  vague  criterion  is  implemented 
in  a  definite  way  in  order  to  permit  illustrative  values  to  be  computed. 

For  a  given  current  threat  state,  z  and  z,  downrange  and  crossrange  distances  {a.b)  to 
the  ship  from  the  threat  are  defined  by 

a  =  jc  COS0  =  -Z‘t/\t\  /c  IX 

h  =  i  sine  =  |xxi|/|i-|  ' 

where  6  is  the  current  look  angle  of  the  ship  as  seen  by  the  threat  (i.e.,  the  angle  between  the 
vectors  -z  and  z).  Recall  that  represents  the  maximum  achievable  lateral  acceleration. 
Then  ibj  =  IxP/a^  is  the  minimum  achievable  turn  radius.  The  domain  of  capability  (described 
in  Chapter  4)  is  the  region 

a  >  0  and  a^+ib-b^)^  >  bf  (5-2) 

in  the  iz,  z)  plane.  If  the  above  condition  is  not  satisfied,  it  is  assumed  that  it  is  not  possible 
to  hit  the  ship.  If  the  condition  is  satisfied,  various  trajectories  to  the  ship  are  possible.  In 
particular,  there  is  a  unique  CTR  trajectory  to  the  ship  with  an  achievable  turning  rate. 

The  LA  trajectories  are  defined  as  follows.  If 

a  >  jbj  and  b  >  bj  (5-3) 

the  trajectory  length  is  taken  to  be  the  larger  of  two  quantities.  The  fu^t  is  the  time  taken  to 
travel  the  distance  a  +  b  minus  the  amount  cut  off  by  rounding  the  comer  with  a  circular  arc  of 
radius  bj .  The  second  is  the  time  taken  to  fly  a  CTR  trajectory  to  the  ship.  There  will  be 
deflnite  regions  of  the  (a,b)  plane  where  one  of  these  alternatives  is  to  be  taken  as  the  LA 
time,  but  these  regions  have  not  been  delineated.  If 

b  <  bi  (5-4) 

the  LA  trajectory  will  consist  of  a  CV  continuation  until  the  latest  opportunity  to  turn,  with 
maximum  achievable  lateral  acceleration,  to  be  able  to  hit  the  ship.  If 

a  <  bj  and  b  >  b^  (5-5) 

the  LA  trajectory  is  the  CTR  trajectory  which  will  hit  the  ship. 
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CHAPTER  6 

COMPARISON  OF  THE  METHODS 


The  prediction  methods  discussed  in  previous  chapters  are  compared  in  Table  6-1.  Each 
column  gives  the  characteristics  of  one  of  the  prediction  methods,  while  each  row  represents  a 
particular  feature  for  all  the  predictors. 


MEASURES  OF  EFFECTIVENESS  (MOE) 


The  MOEs  are  intended  to  quantify  the  ability  of  a  defensive  missile  or  bullet  to  kill  the 
threat.  The  following  parameters  measure  the  effectiveness  of  the  prediction  scheme  in  a  given 
scenario: 

Look  Angle  y(t,tj) 

Heading  Error  ^{t,  t^) 

CPA  Distance  6  ( tj) 

where  is  the  launch  time  of  the  defensive  missile,  and  t  is  time  during  the  engagement.  It 
is  noted  ^t  these  MOEs  are  related  to  a  particular  defensive  missile  to  be  fired  at  the  threat  and 
will  depend  on  the  type  of  trajectory  assumed  by  the  missile.  Thus,  they  are  defmed  only  in  the 
context  of  such  a  defensive-missile  system.  In  the  examples  discussed  below,  the  defensive 
missile  was  a  bullet  flying  a  CV  trajectory.  The  Look  Angle  and  the  Heading  Error  vary  along 
the  missile  trajectory,  and  thus  depend  on  t,  whereas  the  missile-to-threat  distance  at  the  closest 
point  of  approach  (CPA  Distance)  depends  only  on  information  available  at  launch  time;  i.e., 
on  the  thrrat  state  for  any  given  prediction  method. 

The  Look  Angle  is  the  angle  between  the  LOS  from  missile  to  threat  and  the  missile 
velocity.  (In  a  more  detailed  study,  the  Look  Angle  would  be  defined  in  terms  of  the  body  axis 
instead  of  the  missile  velocity.)  The  Look  Angie  is  an  important  consideration  in  the  case  of 
a  missile  with  terminal  homing.  If  the  Look  Angle  is  too  large  at  acquisition  time,  acquisition 
will  be  impossible  because  of  seeker  limitations. 

The  Heading  Error  is  defmed  here  as  the  angle  between  the  velocity  of  the  defensive  missile 
and  the  direction  from  the  current  missile  position  to  the  Perfect-Prediction  PIP  (PPPIP).  While 
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the  PIP  indicates  the  projected  intercept  point  based  on  the  predicted  threat  trajectory,  the  PPPIP 
is  based  on  the  true  threat  trajectory.  It  is  the  threat’s  true  position  at  the  time  when  the  missile, 
having  been  launched  at  time  t  ^ ,  could  just  arrive  provided  it  were  shot  in  the  appropriate 
direction  (which  is  not  necessarily  the  direction  it  was  shot)  as  shown  in  Figure  6-1.  The 
heading  error  is  approximately  the  amount  by  which  the  missile  must  turn,  during  the  time 
remaining  until  intercept,  to  obtain  zero  miss.  The  word  "approximately"  is  used  here  because 
the  PPPIP  is  defined  in  terms  of  where  the  intercept  would  have  occurred  had  the  missile  been 
fired  in  the  appropriate  direction.  Even  if  the  intercept  can  be  achieved,  it  is  at  a  slightly 
different  point  than  PPPIP,  because  the  total  arc  length  from  launch  -.int  to  PPPIP  for  the  actual 
shot  (assuming  that  it  is  capable  of  reaching  the  PPPIP)  is  different  from  the  arc  length  based 
on  the  hypothetical  well-aimed  shot. 

The  Closest-Point  of  Approach  (CPA)  distance  is  the  shortest  distance  between  the  defensive 
missile  and  the  threat  that  will  occur  (based  on  the  actual  missile  and  threat  trajectories). 


TEST  SCENARIO 

A  threat  trajectory  is  considered  that  lies  in  a  single  plane  and  consists  of  several  CV 
segments  joined  by  circular  arcs.  The  threat  speed  is  a  constant  Mach  2  throughout,  and  the 
lateral  acceleration  in  all  turns  is  6g,  as  shown  in  Figure  6-2.  The  various  track  prediction 
schemes  are  illustrated  by  showing  the  predicted  tracks  determined  at  various  points  along  the 
actual  threat  trajectory.  Each  of  these  predicted  tracks  branch  from  the  threat  track  at  a  certain 
launch  time.  Each  predicted  track  ends  at  the  PIP  determined  for  a  hypothetical  flight  of  a 
constant- velocity  bullet  shot  from  the  origin  at  launch  time.  Defense  with  a  higher-speed  bullet 
results  in  PIPs  lying  closer  to  the  point  at  which  each  predicted  track  branches  from  the  actual 
threat  track. 

Figures  6-3  through  6-8  show  the  predicted  CV,  MM,  EA  and  TPN  trajectories  for  various 
parameter  choices.  Note  that  the  scale  varies  somewhat  between  figures.  It  is  seen  that  without 
the  goal-based  hypothesis,  the  predicted  tracks  can  diverge  sharply  from  the  true  threat 
trajectory,  as  the  CV  tracks  of  Figure  6-3.  Figure  6-9  shows  the  PIPs  for  the  CV  predictor  and 
corresponding  actual  synchronous  threat  positions,  joined  by  lines.  Figure  6-10  shows  a 
comparison  of  the  PIPs  for  the  APN,  CV,  and  MM  Predictors,  which  again  illustrates  the 
smaller  divergence  of  the  goal-based  PIPs  from  the  true  threat  trajectory.  This  is  not  surprising 
in  view  of  the  fact  that  this  particular  threat  track  is  goal-based. 

Figure  6-11  shows  the  variation  of  the  MM  acceleration  as  function  of  launch  time  for  two 
values  of  s .  The  sharp  rise  in  the  curve  for  s  =  5kin,  near  the  end  of  the  trajectory,  occurs  at 
the  point  where  the  threat  is  no  longer  able  to  reach  the  ship.  Figures  6-13  and  6-14  show  a 
comparison  of  the  variation  of  the  various  terminal  time  parameters  for  a  pair  of  values  for  the 
limiting  acceleration  and  s .  It  is  seen  that  the  MM  times  are  not  very  different  in  the  two  cases. 
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"y  X  Trajectories 

/  y  a  —  True  Threat  Trojectoiy 
/  y  b  —  Predicted  Threat  Trajectory 
/  c  —  Actuol  Missile  Trajectory 

y  d  —  Missile  Homing  Trojectory 

/  e  —  Original  Guided  Missile  Trajectory 

f  —  Unknown  Best  Missile  Trajectory 

Special  Points 

—  Origin  (Launch  Platform) 

—  PIP 

—  (Unknown)  Perfect  Prediction  PIP  (PPPIP) 

—  Actual  Intercept  Point 

—  Missile  Position  at  Acquisition 

—  Threat  Positian  at  Predictian  Time 


FIGURE  6-1.  TRUE  AND  HYPOTHETICAL  THREAT  AND  MISSILE  TRACKS 
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HGURE  6-2.  HYPOTHETICAL  THREAT  TRAJECTORY  FOR  THE  STUDIES 
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HGURE  6-3.  CV  TRACKS  FOR  A  MACH-2  BULLET 


*10^  MM  Tracks  with  s=5  km 


X  (m) 


xIO* 


HGURE  6-4.  MM  TRACKS  WITH  s=5  km,  MACH-1  BULLET 


6-5 


NSWCDD/TR-92/445 


xIO*  EA  Tracks  for  01=89 
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FIGURE  6-5.  EA  TRACKS  WITH  a^  =  8g,  MACH-1  BULLET 


xIQA  EA  Tracks  for  ol=15g 


HGURE  6-6.  EA  TRACKS  WITH  =  15  g,  MACH-1  BULLET 


6-6 


y  (m) 


NS WCDD /TR -92/445 


6-7 


NSWCDD/TR-92/445 


this  difference  being  illustrated  in  Figure  6-13. 

Figures  6-15  to  6-20  show  the  variation  in  look  angle  and  heading  error  for  some  of  the 
predictors.  It  is  seen  that  at  time  when  the  threat  is  carrying  out  even  a  modest  maneuver  these 
MOEs  assume  very  disadvantageous  (large)  values.  Figures  6-21  and  6-22  compare  the  various 
predictors  according  to  their  variation  in  look  angle.  Figures  6-23  and  6-24  compare  the  same 
predictors’  time  profile  of  heading  error.  Figures  6-25  and  6-26  coixqiare  various  predictors’ 
time  profile  of  miss  distance.  Figure  6-27  illustrates  the  superiority  of  the  (goal-based)  PN 
predictor  over  the  state-based  predictors,  as  is  seen  by  the  lower  magnitudes  of  both  the  look 
angle  and  heading  error  as  the  threat  nears  its  goal.  PIP  tracks  for  two  of  the  PN  predictors 
are  shown  in  Figure  6-28. 


FIGURE  6-9.  PIPs  AND  SYNCHRONOUS  THREAT  POSITIONS  FOR  CV  PREDICTOR 
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xIO*  PIPa  for  o  Mach— 3  bullet 


FIGURE  6-10.  COMPARATIVE  PIP  POSITIONS  (o  =  CV,  +  =  MM,  *  =  APN) 


Minimax  Acceleration,  s««0  (solid),  s»5  km  (dashed) 


FIGURE  6-11.  MM  ACCELERATION  FOR  TWO  VALUES  OF  s 
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MM  time  differance.  b~0  via  o  via  a~5  km 


FIGURE  6-12.  TIME  DIFFERENCE  BETWEEN  THE  MM  ARRIVAL  TIMES  FOR  TWO 
VALUES  OF  s.  (THE  ACTUAL  MM  TIMES  ARE  SHOWN  IN  FIGURES  6-13  AND  6-14) 


Time  Parameters  for  al=15g,  8=5000m 


FIGURE  6-13  VARIATION  OF  TIME  PARAMETERS  WITH  LAUNCH  TIME 
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Tima  Porametera  for  aiB2g,  a“0 


SURE  6-14.  VARIATION  OF  TIME  PARAMETERS  WITH  LAUNCH  TIME 


Look  Anglo  (solid)  ond  Heading  Error  (dashed)  for  CV  Predictor 


nOURE  6-15.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  CV  PREDICTOR 
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Look  Angle  (eolid)  and  Heading  Error  (doehed)  for  CA  Predictor 


HGURE  6-16.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  CA  PREDICTOR 


HGURE  6-17.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  CTR  PREDICTOR 
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Look  Anglo  (aoUd)  ond  Hooding  Error  (doshod)  for  MM  Prodictor 


O  20  40  eO  80  100  120  140  160  180  200 
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FIGURE  6-18.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  MM  PREDICTOR 

Look  Anglo  (lolid)  ond  Heading  Error  (dashed)  for  EA  Predictor 
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nOURE  6-19.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  EA  PREDICTOR 
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Look  Angle  (solid)  ond  Heading  Error  (dashed)  for  ED  Predictor 


HGURE  6-20.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  EDTR  PREDICTOR 


Look  Angle  for  CV  (solid),  CA  (doshed),  and  CT  (dotted)  Predictors 


HGURE  6-21.  LOOK  ANGLE  FOR  CV,  CA,  AND  CTR  PREDICTORS 
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Look  Angle  for  ED  (solid),  MM  (dashed),  and  EA  (dotted)  Predictors 


HGURE  6-22.  LOOK  ANGLE  FOR  EDTR,  MM,  AND  EA  PREDICTORS 


Heading  Error  for  CV  (solid),  CA  (dashed),  ond  CT  (dotted)  Pred. 


FIGURE  6-23.  HEADING  ERROR  FOR  EDTR,  MM,  AND  EA  PREDICTORS 
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Hsoding  Error  for  EO  (solid).  MM  (doshad),  and  EA  (dotted)  Pred. 


IGURE  6-24. 


HEADING  ERROR  FOR  CV,  CA,  AND  CTR  PREDICTORS 


X10-*  CV  (solid),  CA  (dashed).  CT  (dotted) 


FIGURE  6-25.  MISS  DISTANCE  FOR  CV,  CA,  AND  CTR  PREDICTORS 
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xIO*  EO  (solid),  MM  (doshsd),  EA  (dotted) 


3URE  6-26.  MISS  DISTANCE  FOR  EDTR,  MM,  AND  EA  PREDICTORS 


nOURE  6-27.  LOOK  ANGLE  AND  HEADING  ERROR  FOR  PN  PREDICTOR 


6-17 


NSWCDD/TR- 92/445 


I 

I 


6-18 


NSWCDD/TR-92/445 


CHAPTER? 

SEA-SURFACE  CONSTRAINED  THREAT  TRAJECTORY  PREDICTION 


The  previous  discussion  has  assumed  that  the  predicted  threat  trajectories  do  not  intersect 
the  sea  surface.  This  is  not  always  the  case.  Diving  threats,  in  particular,  will  tend  to  generate 
predicted  trajectories  intersecting  the  sea  surface.  Modifications  are  proposed  here  so  that  the 
predicted  trajectories  will  lie  entirely  above  the  sea  surface.  Similar  modifications  are  needed 
to  redefine  some  of  the  Terminal  Parameters  so  that  they  be  based  on  "dry"  trajectories.  Only 
the  planar  trajectories  are  studied  in  this  chapter,  and  for  these  the  sea  is  to  be  represented  by 
the  line  where  the  sea  surface  plane  intersects  the  flight  plane.  This  will  be  referred  to  as  the 
"sea  line." 

Let  zi  be  a  unit  vector  in  the  flight  plane,  perpendicular  to  the  sea  line  and  directed  so  that 
0.  Because  the  threat  will  never  cross  the  sea  line,  a  remains  constant  and  can  be 
determined  from  initial  conditions. 


MM  TRAJECTORY 


The  constraint  is  adjusted  if  the  end  of  the  circular-arc  segment  is  below  the  sea  line.  This 
occurs  if 


sin0  . 


0) 


^0^ 


1 


-COS0 c  \ 
0)2 


•a  <  0 


(7-1) 


With  the  definitions  of  Chapter  4,  the  vector  quantity  within  the  parentheses  is  clearly  the  threat 
position  at  the  transition  point  between  the  CTR  and  CV  segments.  If  the  situation  (7-1)  arises, 
ideally,  a  trajectory  avoiding  the  sea  line  but  reaching  the  origin  with  minimal  lateral- 
acceleration  magnitude  should  be  found.  The  following  simplifrcation  is  hypothesized  to  provide 
such  a  trajectory,  but  no  demonstration  of  this  has  been  found.  The  idealized  trajectory  consists 
of  one  circular  arc  in  the  flight  plane  with  CV  segments  before  and  after  it.  The  trajectory 
having  these  properties,  plus  the  minimum  possible  lateral-acceleration  magnitude,  is  the  one  that 
is  determined  by  the  following  scheme.  Define  the  "dive  point, " 


- rJ 

fi'X 


(7-2) 
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where  the  threat  would  impact  the  sea  line  if  a  CV  trajectory  were  taken.  Then,  at  a  distance 

a  =  min[lr^| -s,  Ix-Xjl]  (7-3) 

before  reaching  the  dive  point;  that  is,  when  the  threat  is  at  the  position 


a  circular-arc  segment  is  begun.  The  total  arc  length  6  in  this  segment  is  given  by 


where  the  radius  of  curvature  is 


COS0'  =  - 


(7-5) 


h  = 


a 

tan  {6'/2) 


(7-6) 


If  h  turns  out  to  be  too  small  (requiring  an  acceleration  magnitude  that  is  too  large),  no  MM 
trajectory  is  possible.  The  center  of  curv’ature  is  at  the  point 


where 


v/a2+jb2  +  ^2 - 2r 


(7-7) 


(7-8) 


is  the  point  on  the  sea  line  where  the  circular-arc  segment  ends  and  the  final  CV  sea-skimming 
run  over  the  distance  s  is  begun.  The  initial  acceleration  in  the  circular  arc  is  given  by 


r  = 


\i\^ 


(7-9) 


Figures  7-1  though  7-5,  show  on  an  arbitrary  scale  some  typical  MM  tracks  with  a  run-in 
distance  of  s  =  2000/n  and  sea  constraint.  In  Figure  7-1,  the  condition  (7-1)  does  not  arise, 
and  the  trajectory  is  not  affected  by  the  presence  of  the  sea.  In  Figures  7-2  through  7-5,  there 
is  a  final  CV  segment,  at  least  as  long  as  s ,  along  the  sea  line,  with  successively  greater  lateral- 
acceleration  magnitudes  required  in  the  CTR  segment. 


Since  no  demonstration  that  the  foregoing  trajectory  does  indeed  minimize  the  maximum 
lateral-acceleration  magnitude  has  been  given,  there  remains  a  possibility  that  a  double-arc 
trajectory,  first  turning  away  from  and  then  toward  the  sea  line,  would  satisfy  the  constraints 
and  result  in  a  lower  maximum  lateral-acceleration  magnitude.  The  geometry  is  a  bit  messy, 
and  this  possibility  has  not  been  investigated. 
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FIGURE  7-2.  EXAMPLE  OF  TRAJECTORY  IN  ir.i)  PLANE 
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MM  Predicted  Track  Constrained  by  Sea  Line  (Dashed  Curve) 
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FIGURE  7-5.  EXAMPLE  OF  TRAJECTORY  IN  (x,i:)  PLANE 


EA  TRAJECTORY 

The  effect  of  the  sea  line  on  the  EA  predictor  is  quite  simple.  If  the  unconstrained  EA 
trajectory  crosses  the  line,  there  is  simply  no  constrained  EA  trajectory  at  all.  This  is  because 
the  limiting  lateral-acceleration  magnitude,  ,  is  specified.  A  sharper  turn  to  avoid  the  sea  line 
is  impossible.  A  simple  check  is  made,  as  in  the  MM  case  discussed  in  the  previous  section, 
but  with  the  computations  based  on  the  p  computed  as  in  Chapter  4. 


PN  TRAJECTORY 

The  predicted  trajectories  can  be  kept  away  from  the  sea  line  by  an  additional  term  in  the 
required  acceleration.  Let  be  a  vector  whose  magnitude  is  the  distance  between  the  threat  and 
the  sea  line  and  in  the  direction  of  b-  This  is  equal  to 

b  =  ia'r)B  (7-10) 

If  the  threat  were  to  fly  CV,  it  would  hit  the  sea  line  in  the  time 
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(7-11) 


It  is  desirable  to  introduce  an  incremental  threat  acceleration  only  if  the  threat  would  otherwise 
cross  the  sea  line.  Let  this  incremental  acceleration  be  denoted  by  t  ■  This  will  be  taken 
perpendicular  to  the  threat  velocity,  and  with  magnimde  that  increases  inversely  with  distance 
from  the  sea  line,  and  directly  with  velocity  component  toward  the  sea  line.  These  requirements 
are  met  by  the  expression 


lij  =  4^  if  t,>  0 

llj  =  0  if 


(7-12) 


The  incremental  threat  acceleration  must  have  the  form 


I 


S 


(7-13) 


in  order  to  lie  in  the  (z;,x)  plane  and  be  perpendicular  to  the  threat  velocity.  The  scalar 
quantity  p  is  determined  in  order  that  the  incremental  acceleration  have  the  appropriate 
magnitude  (7-12).  Accordingly, 


V-  = 


{U'Z) 


(r-±) 


1 

2 


(7-14) 


However,  improvement  is  still  needed.  The  above-described  incremental  acceleration  can 
be  appreciable  even  in  cases  where  the  threat  would  not  cross  the  sea  line,  especially  in  cases 
where  the  threat  velocity  is  directed  nearly  normal  to  the  sea  line.  (It  does  not  matter  if  the  sea 
line  is  crossed  at  the  point  occupied  by  the  ship.)  Therefore,  it  is  advantageous  to  multiply 
by  a  factor  that  is  relatively  small  for  normal  incidence  (of  the  threat  velocity  toward  the  sea 
line)  and  increases  to  unity  for  glancing  incidence.  There  is  considerable  arbitrariness  in 
specifying  such  a  factor,  but  here  it  has  been  found  convenient  to  use  sin^6 ,  where  5  is  the 
angle  between  the  threat  velocity  and  xj.  This  angle  is  given  by  the  relation 


Thus, 


cos  6 


(7-15) 


sin®6  = 


\i\^ 


(7-16) 


This  factor  reduces  unwanted  modification  of  the  PN  law  in  cases  where  no  modification  is 
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needed.  In  addition,  there  is  no  a  priori  justification  for  applying  the  incremental  acceleration, 
as  just  defined,  with  unit  gain.  Let  v  be  an  arbitrary  gain,  to  be  adjusted  to  achieve  the  desired 
effect.  The  incremental  acceleration  is  then  represented  by 


1  - 

(a-x)^V 

‘L 

_  1 
'  2 

nz 

JL 

■r" 

J 

(7-17) 


This  expression  may  appear  a  bit  complicated,  but  it  serves  the  purpose  and  is  easy  to  compute. 
Figures  7-6  to  7-9  show  the  resulting  APN  and  MPN  trajectories  for  various  values  of  v .  The 
value  V  =  0  results  in  the  unmodified  PN  law  and  in  these  cases  results  in  a  trajectory  that 
crosses  the  sea  line.  Small  positive  values  (on  (0,.3))  cause  numerical  difficulties  near  the  sea 
line  and  are  not  shown.  The  results  indicate  that  the  value  v  =  .  5  is  appropriate.  Figures  7-6 
and  7-7  show  two  examples  of  Adler  PN,  and  Figures  7-8  and  7-9  show  two  examples  of  MPN. 
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FIGURE  7-6.  EXAMPLE  OF  TRAJECTORY  IN  (x,i:)  PLANE 
(ARBITRARY  SCALE) 
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CHAPTERS 

CONCLUSIONS 


The  state-based  predictors  are  easily  deceived  by  maneuvering  threats  into  producing  PIPs 
that  lie  far  from  the  actual  threat  track  whose  destination  is  the  ship.  These  state-based 
predictors  are  distinguished  by  the  degree  of  confidence  placed  on  the  estimated  threat 
acceleration  and  jerk.  With  minimal  confidence  in  the  estimates  of  acceleration  and  jerk,  one 
can  do  no  better  than  a  CV  linear  extrapolation  of  the  initial  threat  position.  Maximal 
confidence  in  these  estimates  leads  one  to  the  Helical  predictor.  Between  these  extremes  are  the 
CTR  and  EDTR  predictors,  which  ignore  the  initial  jerk,  with  the  former  applying  the  initial 
lateral  acceleration  throughout  the  prediction  while  the  latter  has  the  acceleration  decay 
exponentially  to  zero. 

The  goal-based  predictors  were  found  to  be  useful  when  the  destination  of  the  threat  is 
somewhat  unambiguous.  All  ignore  the  initial  estimated  threat  acceleration  except  the  TPN, 
which  may  be  used  when  this  value  is  deemed  reliable  and  important.  The  large  prediction 
errors  typically  generated  by  the  state-based  predictors  are  considerably  reduced.  However,  the 
state-based  predictors  are  usually  more  accurate  for  short  times. 

Most  of  the  predictors  used  the  constant-speed  constraint.  While  this  constraint  may  not  be 
realized  by  all  threats,  taking  speed  variation  into  account  requires  threat  models  that  are 
considerably  more  complicated.  The  only  exceptions  were  the  CA  and  the  CLLA  predictors, 
in  which  the  speed  can  vary.  The  CA  predictor  is  not  recommended  except  for  short-term 
predictions,  llie  CLLA  predictor  may  be  useful  if  the  variation  in  speed  can  be  estimated. 
Methods  for  adjusting  the  initial  threat  state  estimates  to  satisfy  the  constant-speed  constraint 
were  outlined. 

A  number  of  Terminal  Parameters  were  defined  in  Chapter  5  for  the  purpose  of  facilitating 
decision  making  by  providing  the  user  with  more  innmediately-usable  information.  For  example, 
the  parameter  the  earliest  arrival  time,  might  be  monitored  for  each  threat  to  ensure  that 
the  target  is  engaged  in  a  timely  manner.  The  other  Terminal  Parameters  similarly  represent 
the  voluminous  essential  information  into  a  more  concisely  manageable  form. 

The  concepts  discussed  in  this  report  need  to  be  developed  further.  The  next  advance  in 
methodologies  for  dealing  with  maneuvering  threats  may  involve  specifics  of  the  different  threat 
types,  as  knowledge  of  their  projected  capabilities  becomes  known.  Future  prediction  schemes 
may  use  target  ID  when  available,  or  characterize  its  behavior  in  an  adaptive  (on  line)  algorithm. 
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